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Abstract 

The muon anomalous magnetic dipole moment (MDM) is calculated in the 
framework of the minimal supersymmetric standard model (MSSM). In this paper, 
we discuss how the muon MDM depends on the parameters in MSSM in detail. We 
show that the contribution of the superparticle-loop becomes significant especially 
when tan/3 is large. Numerically, it becomes O(10 -8 — 10 -9 ) in a wide parameter 
space, which is within the reach of the new Brookhaven E821 experiment. 
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1 Introduction 



Supersymmetry (SUSY) [I]] is one of the most attractive candidates of the new physics 
beyond the standard model. In SUSY models, quadratic divergences are automatically 
canceled out, and hence SUSY may be regarded as a solution to the naturalness prob- 
lem 0. In addition, precision measurements of the gauge coupling constants strongly sug- 
gest SUSY grand unified theory (GUT) [||]. Contrary to our theoretical interests, however, 
evidences of SUSY (especially, superpartners) have not discovered yet, and hence super- 
partners are fascinating targets of the forthcoming high energy experiments like LEP II, 
LHC and NLC. 

Even if we do not have high energy colliders, we can constrain SUSY models by 
using precision measurements in low energy experiments. This is because superparticles 
contribute to low energy physics through radiative corrections. Especially, superparticles 
are assumed to have masses of the order of the electroweak scale, and hence their loop 
effects may become comparable to those of W ± - or Z-boson propagations. Therefore, low 
energy precision experiments are also very useful to obtain constraints on SUSY models. 

One of the quantities which are measured in a great accuracy is the muon anomalous 
magnetic dipole moment (MDM), = |(g — 2) M . At present, the muon MDM is measured 
to be U 

a™? = 1165923(8.4) x 1(T 9 . (1) 

On the other hand, the standard model prediction on a M is given by H 

a s ™ = 116591802(153) x KT 11 , (2) 

which is completely consistent with experimental value. (For a review of the calculation 
of af 1 , see also Ref. §.) 

Because of the great accuracy of a^ xp and a^ M given above, we can derive a constraint 
on SUSY models from the muon MDM. Furthermore, the new Brookhaven E821 exper- 
iment [0] is supposed to reduce the error of the experimental value of to 0.4 x 10~ 9 , 
which is smaller than the present one by a factor ~ 20. The accuracy of the Brookhaven 
E821 experiment is of the order of the contribution of the W - and Z-boson loop, which 
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means we may have a chance to measure the SUSY contribution to the muon MDM by 
that experiment. 

In fact, there are several works in which the muon MDM is calculated in the frame- 
work of SUSY models M, |S|, [TO] . Especially, Chattopadhyay and Nath recently pointed 



out that the muon MDM is a powerful probe of the models based on supergravity if 
tan/? is large However, most of the recent works assume the boundary conditions 

on the SUSY breaking parameters based on the minimal supergravity, and/or radiative 
electroweak symmetry breaking scenario, and hence it is quite unclear for us how the 
SUSY contributions to the muon MDM depend on the parameters in MSSM. Thus, the 
aim of this paper is to clarify it, and to investigate the behavior of the muon MDM in 
the framework of MSSM. The mass matrices and mixing angles among the superparticles 
have model dependence even if we assume the boundary condition based on the minimal 
supergravity, and hence we believe that it is important to analyze the muon MDM in a 
more general framework of the SUSY standard model. 

In this paper, we investigate the SUSY contribution to the muon MDM in the frame- 
work of MSSM as a low energy effective theory of SUSY GUT [PH| . The organization 
of this paper is as follows. In the next section, we introduce a model we consider. In 
Section || we show analytic forms of the SUSY contribution to the muon MDM, Aa^ USY . 
In Section f|, typical behavior of Aa^ USY is discussed. In Section [5], some numerical results 
are shown. Section |6| is devoted to discussion. 

2 Model 

First of all, we would like to introduce a model we consider, i.e. MSSM as a low energy 
effective theory of SUSY GUT. All the fields we use in our analysis are 

l L = fa u), ii% H x = (#r Hi), H 2 (3) 
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where II (2*, — ~) and \i c R (1, 1) are left- and right-handed leptons in the second generation, 
while two Higgs doublets are represented as H\ (2*, — |) and H 2 (2, |). (We denote the 
quantum numbers for the SU(2) L x U(l)y gauge group in the parenthesis.) The Higgs 
doublets H\ and H 2 are responsible for the electroweak symmetry breaking, and hence 
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their vacuum-expectation values are constrained as (Hi) 2 + (H2) 2 — (174GeV) 2 in order 
to give a correct value of the Fermi constant. On the other hand, the ratio of the vacuum- 
expectation values of two Higgs doublets is a free parameter in MSSM, which we define 
t&n(3={H 2 )/{H x ). 

Relevant part of the superpotential of MSSM is given by 

W = y^n c R lL„H w + fi H H ha H«, (4) 

where y^ is the Yukawa coupling constant of muon, \xn the SUSY invariant Higgs mass 
and e a/3 the anti-symmetric tensor with e 12 = 1. Using the superpotential given above, 
F-term contribution to the lagrangian is obtained as 

C F = - J d 2 6W + h.c. (5) 

Furthermore, soft SUSY breaking terms are given by 

£ so ft = -m 2 J*J L - m 2 R jl c R jl R - {A ll e ap pL R l L ^H2,p + h.c.) 

--(m G2 WW + m G1 BB + h.c). (6) 

Here, l^, jl R , W and B represent left- and right-handed sleptons in second generation, 
and gauginos for SU(2)l and U(l)y gauge group, respectively. Gaugino masses ma and 
are related by the GUT relation; 

777 G2 = 3 777G1 

where g\ and gi are the gauge coupling constant of SU(2)l and U(1)y gauge group, 
respectively^ 

Here, we should comment on a flavor mixing in slepton mass matrices. If there is large 
flavor mixings in the slepton mass matrices, all the sleptons contribute to the muon MDM. 
However, flavor mixing in the slepton mass matrices may be dangerous, since it induces 
lepton flavor violation processes such as \i — > ej, r — > /ry and so on. Especially, the 
mixing among first and second generation is severely constrained from \x — > especially 



when tan/3 is large ||13||. On the other hand, the constraint on the mixing of the second 



^^The GUT relation given in eq.(0) holds in general if the gauge groups are unified in a larger group |H 
Therefore, we are not depending on specific model of GUT. 
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and third generation is not so stringent. In this paper, for simplicity, we assume that the 
flavor mixing in the slepton mass matrix is not so large, and that it does not affect the 
following arguments. A comment on the case with the flavor mixing is given in Section |6|. 

Once we have the MSSM lagrangian, we can obtain mass eigenvalues and mixing 
matrices of the superparticles. The mass matrix for the smuon field is given by 



where 



m 



LR 



Mi 




LR 



m 
m 



LR 
2 



(8) 



m l + m l cos 2/? ( sin 2 9w 



m\ — mf cos 2(3 sin 2 9w, 



(9) 
(10) 

(11) 



The mass matrix M~ can be diagonalized by using an unitary matrix Up as 

(UtMp^AB = m\ A b AB (A, 5 = 1,2), 



(12) 



where m^A is the mass eigenvalue of the smuon. Notice that, in our case, off-diagonal 
element of the mass matrix given in eq.(^|) is substantially smaller than the diagonal 
elements, and hence and mp,R almost correspond to the mass eigenvalues. The mass 
of the sneutrino, mo, is also easily obtained as 



m,- 



m L H — m| cos 2/3. 



(13) 



Next, we derive the mass matrices for neutralinos and charginos. For neutralinos, the 
mass terms are given by 



— (iB iW 3 H° # 2 °) 
/ 



-men. 




+h.c, 





-m G2 

1 



n/2 



\l H 

Mh / 



/ iB \ 
iW 3 
Hi 
\ H° 2 J 



(14) 
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where H\ and H2 represent the higgsino field. Then, we can find an unitary matrix U x o 



which diagonalize the mass matrix given above. Denoting the mass matrix given in eq.(14) 
as M x o, mass eigenvalues m x o x for the neutralino Xx * s gi yen by 

{UloM x oU x o)xY = m x0x 8 XY (X,F = l-4), (15) 

where m x o x < m x Q Y if X < Y. 

Similarly, mass terms for the charginos are given by 

with W = —^{W 1 =F iW 2 ). The mass matrix given in eq.(plj), which we denote M x ±, 
can be diagonalized by using two unitary matrices, U x + and U x - ; 

{Ul + M x ±U x -) XY = m x ± x 8 XY (X,Y = 1,2), (17) 

where m x ± x represents the mass eigenvalue of the chargino field. 

With the coupling constants and mixing matrices given above, we can write down 
muon-neutralino-smuon and muon-chargino-sneutrino vertices. Denoting the mass eigen- 
states of the smuon, neutralino and chargino as Jxa, Xx an d Xxi respectively, the interac- 
tion terms are given by 

C mt = £ fi(N% x P L + N* x P R )x x fiA + E HCxPl + C*Pn)xtv + h.c. , (18) 

AX X 

where P L = |(1 - 75), Pr = |(1 + 7s) and 

N L AX = -V^ohxiU^LA-^iU^xiU^RA, (19) 

Kx = -^(^x03x(^)^ + ^2(^ x o) 2X (^) LA + ^ 1 (f/ x o) lx (f/ A ) LA , (20) 

c L x = y»{u x -) 2X , (21) 

C§ = -92(U X+ ) 1X . (22) 

By using the interaction terms given in eq. (|T8|) , we calculate the SUSY contribution to 
the muon MDM. 
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3 Analytic formulae 

Now, we are in position to calculate the SUSY contribution to the muon MDM. What we 
have to calculate is the "magnetic moment type" operator, which is given by 

£mdm = ^F 2 fia p xfiF p \ (23) 

Here, e is the electric charge, m p the muon mass, a p \ = f [7 p ,7a], F p \ the field strength 
of the photon field and F 2 the magnetic form factor. The muon anomalous magnetic 



moment, a p , is related to F 2 as 



a p = F 2 . (24) 



Thus, by calculating magnetic form factor in the framework of MSSM, we can have SUSY 
contribution to the muon MDM. 

In SUSY model, there are essentially two types of diagrams which contribute to a p , 
i.e. one is the neutralino (x°)-smuon [fx) loop diagram (Fig. |]a) and the other is the 
chargino (x )-sneutrino (p) loop diagram (Fig. |T]b); 

Aa susv = Aa x°A + Aa x^ (25) 
Here, contribution from the x°~fi diagram, Aa^°^, is 

A< A = ^E{" m ,( N AX^ x + N^ x N^ x )ml A .Um 2 xQx , m| A , m| A , mJ A , m\ A ) 

AX 



+m x0x N^ x N^ x J 4 (m 2 x o x ,m 2 x o x , m 2 - A , m| A )} 



1^1 m,. 



m y —^t r-AN L AX ^ L AX + N% X N* X 

ixX Qmi A (i- XAX y y AX AX AX AX ' 

2 i o™3 aJi 



x (1 - 6x A x + 3x AX + 2x AX - Qx AX In x A x) 

lN L AX N^ x {\ - x\ x + 2x AX hiXAx)), (26) 



m x°x f l AT R n ^ 
m lA{ l ~ x Axf 



where we are using mass eigenstate basis of x° an d (and that of \ in deriving eq. 
Here, x AX = m x o x / m2 i A , and we define the functions In and Jjy as 

I N {m 2 ,...,m 2 N ) = / (2 t^- (fe2 _ mf) . 1 . (F _ <r (27) 

J N (m 1 ,---,m N ) = / /0 \a- 7T2 2\ 772 2V ( 2 °) 
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Some useful formulae concerning the functions In and Jn are shown in Appendix |A]. 
Contribution from the x & loop diagram is also easily calculated, and the result is given 
by 



where xx 



in 



£ [m,{C L x C L x + C7|C|){ J 4 K ±X , m x±x , m x±x , m x±x ) 

i 2 t / 2 2 2 2 2\ T /2 2 2 2\1 

+m i> J 5 (m x ± x , m x±x , m x±x , m x±x , m p ) - J 4 (m x±x , m x±x , m x±x , m^)} 
-2mxC^C| J 4 (m x±x , m x±x , m x±x , m~) 



16vr 2 

x 



^\3m?(l 



1 + -xx - 3x x + + 3x x In 



3m 



m|(l -x x ) 3 



4 1 2 2 , 

+ -x x + -lnxx 



(29) 



m: 



/m|. 



4 Behavior of the SUSY contribution to the muon 
MDM 

Before evaluating the SUSY contribution to the muon MDM numerically, we would like 
to discuss the behavior of Aa^ USY , especially in large tan/? case. As we will soon see, 
|Aa^ USY | becomes large as tan/3 increases. Thus, the discussion about the large tan/? 
case will be helpful for us to understand the behavior of Aa^ USY . 

For this purpose, it is more convenient to use the mass insertion method to calculate the 
penguin diagrams rather than working in the mass eigenstate basis of the superparticles 
which is used in the previous section. In the case where tan/5 is large, five diagrams 
dominantly contribute to Aa^ USY , which are shown in Fig. |[ Their contributions are 
given by 

= gl m l m GiHH tan 13 
AaJJ 2 = -gfmlm G1 fj, H tan/3 

X {M m GV m Gl^H^ m lw m lR) + M m ll^H^H, m lR, m l R )}, (31) 
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Aajf = -glmlm G iiiH tan (3 

X {M m Gl, m Gl^H, m lL, m lL) + M m Gl^H^H, m lL, m lL)}i (32) 

AaJ^ 4 = -^g 2 mlm G2 fi H tan f3 

x{j5{m 2 G2 ,m 2 G2 ,fi 2 H ,ml L ,ml L ) + J 5 {m 2 G2 , fi 2 H , fi 2 H ,ml L ,ml L )}, (33) 
Aa^ = glm 2 jn G 2^Htaji(3 

„f n r / 2 2 2 2\ 7/2 2 2 2 2\ 

x |2/ 4 (m G2 , m G2 , /i^, m f ) - J 5 {m G2 , m G2 , fi H , m P ) 

+2I 4 (m 2 G2 , ii 2 h , H 2 H ,ml) - J 5 (m 2 G2j m 2 , m~)}. (34) 

Here, eqs.(^) - ( |33| ) are y 0- /^ 1°°P contributions, while eq.(|34"D represents the x^-v loop 
one. By using these expressions, the SUSY contribution to the muon MDM is approxi- 
mately given by 

Aap ~ Aa* 1 + Aaf + Aaf + A< 4 , (35) 
Aaf > ~ AaJ. (36) 

Notice that the SUSY contribution to the muon MDM given in eqs. (|30D - (0) approx- 
imately correspond to the terms which are proportional to N L N R or C L C R {i.e. terms 
which have a chirality flip in internal fermion line) in the exact formulae given in eqs. (|26"D 
and (H). 

The first thing we can learn from the above expressions is that all the terms given in 



eqs. (p0|) - (|34"D are proportional to tan/3 P, |i~0f. This is due to the fact that the chirality 



is flipped not by hitting the mass of the external muon but by directly hitting the Yukawa 



coupling. This mechanism also occurs in the case of the lepton flavor violations fL3 |. 
Thus, |Aa^ USY | becomes large as tan/3 increases, and we obtain severer constraint on the 
parameter space as tan f3 gets larger. 



The second point we should mention is that the relation between the sign of Aa 



SUSY 



and those of parameters in MSSM. The dominant SUSY contribution given in eqs.( 
( |3~3| ) are all proportional to m G n H tan j3 (with m G = rn G i,m G2 being gaugino mass). Thus, 
if we change the sign of this combination, Aa^ u also changes its sign. Furthermore, in 
the case where we assume GUT relation on the gaugino masses, we checked that Aa^ 1 
or Aa^ dominates over other terms (Aa^ 2 , Aa^ 3 , Aa^ 4 ) in most of the parameter space. 



S 



Here, both Aa^ 1 and Aa^ have the same sign as the combination ma^H tan (3. Therefore, 
Aa^ USY becomes positive (negative) when the sign of the combination vtiq^h tan (3 is 
positive (negative). In the next section, we will see that this relation really holds as a 
result of numerical calculations. 

Furthermore, we comment here that the contribution of x^-v loop diagram dominates 
over that of the loop ones if all the masses of the superparticles are almost degenerate. 
For example, let us consider the extreme case where all the masses for the superparticles 
(max, mc2, Hh, n^fiL, n^fiR, mp) are the same. Denoting the masses of the superparticles 
?^susy, contributions of the x°~fi an d X^-v loop diagrams to the muon MDM is given by 

Aa^ * Aaf + Aaf + Aaf + A< 4 

'" {gl-gl) tan/3, (37) 



192tt 2 m| USY 



Aaf° ~ AaJ 



* ^tan/3. (3* 



327T 2 mg USY 



From the above expressions, we can see that the x v loop contribution is substantially 
larger than that of x°-f i loop. Thus, x^-v loop gives dominant contribution, as in the 



case of minimal SUSY GUT based on the minimal supergravity jnj . However, we should 
note here that the x ± '^ 1°°P dominance does not hold in general. In the next section, 
we will see the SUSY contribution to gnificantly depends on the right-handed 

smuon mass m^R in certain parameter regions. 

5 Numerical Results 

In this section, we numerically estimate Aa^ USY by using eqs.(^) and (|29|). As we men- 
tioned before, there are essentially six parameters on which Aa^ USY depends, i.e. SU(2) 
gaugino mass mc2,0 left- and right-handed smuon masses ra^L and m^R (which essen- 

2 If tug or hh is small, this relation does not hold. This is mainly because that the mass insertion 
method breaks down in such a case. Furthermore, in such a case, we cannot ignore Aa^ 2 or terms which 
is not proportional to tan/3 (i.e. terms which are proportional to N L N L , N R N R , C L C L and C R C R in 



the exact formula given in eqs.(|26|) and (|29|)). In that case, the sign of mcfJ-H tan (3 is not directly related 
to that of Aa^ USY . 

3 Gaugino mass for U(1)y gauge group is determined by the GUT relation (0). 
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tially correspond to the soft SUSY breaking parameters m 2 L and m R , respectively), SUSY 
invariant Higgs mass fijj, ratio of the vacuum-expectation values of the two Higgs dou- 
blets, tan/5 = (H2) / (Hi), and SUSY breaking A-parameter for the smuon, A M . However, 
especially in the large tan/3 region where SUSY contribution to Aa^ USY may become 
significantly large, Aa^ USY is not sensitive to if ~ 0(^/1^)- This is because A^ al- 
ways appears in expressions in the combination of (A^ + y^/Xfftan/?), as shown in eq.(|ll]). 
Therefore, in our analysis, we take A^ = 0.[] Then, we take the other five parameters as 
free parameters and calculate the SUSY contribution to the muon MDM for a given set 
of parameters. 

First, we show the SUSY contribution to the muon MDM for fixed values of m^R 
and nifiL in the [iH-mG2 plane. In Fig. |3], we plotted the results for rrijXR = lOOGeV 
and tan/3 = 30. Here, the left-handed smuon mass is taken to be lOOGeV (Fig. |a), 
300GeV (Fig. |b) and 500GeV (Fig. |c). The results for the cases of m AiJ = 300GeV and 
m/m = ITeV are also shown in Fig. ^ and Fig. |^, respectively. As we can see, if we take a 
smaller value of m^R, the SUSY contribution to the muon MDM is enhanced in the large 
fin region. This can be easily understood if we think of the fact that Aa^ 1 gives a large 
contribution in such a parameter region. 

Furthermore, by choosing the right-handed smuon mass m^R so that |Aa^ USY | is min- 
imized, we obtain the lowerbound on the SUSY contribution to the muon MDM. The 
results are shown in Fig. [| Here, we assume 45GeV < m^R < ITeV. The lowerbound 
is obtained from the negative search for the smuon while the upperbound is due to 
the naturalness point of view. In fact, the results are insensitive to the upperbound if we 
take the upperbound larger than about ITeV, since the effects of the right-handed smuon 
decouple when we take m^R — ► 00. 

Here, we would like to discuss the behavior of Aa^ USY shown in Fig. [| As can 

be seen, Aa^ USY changes its behavior at around ~ m^i,- This can be understood 

in a following way. In the case of \hh\ ~ m-fiL, AaJ^ 1 and Aa^ 2 almost cancels out 

and Aa^ USY becomes insensitive to m^R. In the case of ^ m^L, m^-dependence of 

4 The supergravity model suggests ~ 0{y lL mf L ) Furthermore, it was pointed out that some 

unwanted minimum appears in the potential of the smuon when | | > O ( 1 ) x y^m^ Jl5[ , which may cause 
cosmological difficulties. We checked that the results are almost unchanged even if we take A^ = ?>y^m,jxL- 
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Aa^ t USY is almost determined by that of Aa^ 1 . Then, |Aa^ USY | becomes smaller as rrijXR 
becomes larger. On the other hand, if \fijj\ ^ m^i, AaJ^ 2 determines the m^R- dependence 
of Aa^ USY . The important point is that the sign of Aa^ 2 is opposite to that of Aa^ 



which gives the dominant contribution. Thus, |Aa^ | gets smaller as m^R decreases. 
In summary, in the case of \[ih\ ~ "Vl, |Aa^ USY | increases as m^R decreases, while in the 
case of | /iff | ^m^L, |Aa^ USY | decreases as m^R gets smaller. 

Notice that some regions of the ^H~ m G2 plane are excluded by the negative search for 
signals of neutralinos or charginos [16, [17|]. In Fig. |7], we show the excluded region for 



tan/3 = 30, i.e. for large tan/3 case.f] Here, we adopt the following constraints fI7fl 



Ar z < 23.1MeV, (39) 

Ar inv < 8.4MeV, (40) 

Br(Z^xlxl) < 5xHT 5 , (41) 

Br(Z^x°2X°2) < 5xl0- 5 , (42) 

where ATz is the partial width of Z-boson decaying into charginos or neutralinos, while 
Ar inv = Tz{Z — > XiXi) represents the neutralino contribution to the invisible width. For 
the constraint on the chargino mass, we consider several cases where the lowerbound on 
the chargino mass is given by 45GeV (LEP), 90GeV (LEP II), and 250GeV (NLC with 
i/i = 500GeV). Comparing Fig. [7| with Fig. |3| - Fig. |6], we can see that the muon MDM 
has a better sensitivity to MSSM than colliders in some parameter space. 

Remember that the SUSY contribution to the muon MDM is approximately propor- 
tional to tan (3. Therefore, even for the case of tan f3 ^ 30, we can read off the approximate 
value of Aa^ USY from Fig. [3] Fig. || For example, the contours for Aa^ USY = 2 x 10~ 9 
in these figures correspond to Aa^ USY ~ 4 x 10~ 9 for the case of tan/? = 60. 

If the new Brookhaven E821 experiment measures the muon MDM with the accu- 
racy of their proposal, it will give a great impact on MSSM. In a large parameter space, 
|Aa^ USY | becomes O(10 -9 ), which is within the reach of the new Brookhaven E821 experi- 
ment. Furthermore, the theoretical uncertainty, which is almost originate to the hadronic 

5 If tan/3 is fairly large (tan/3^5), mass matrices of the charginos and neutralinos become almost 
independent of tan/3. In this case, the constraint is insensitive to tan /?. We would like to note here that, 
if tan (3 is not so large, in our convention, the constraint becomes severer for the case of hhitt,G2 > 
rather than < 0, as shown in Refs. Cq, O. 
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uncertainty, is also expected to be decreased due to better measurements of the cross 
section of e + + e~ — > hadrons at low energies. Thus, the muon MDM should be regarded 
as a good probe of MSSM. In particular, the Brookhaven E821 may be able to see the 
signal of MSSM even in the case where we cannot find any superparticle by NLC with 
yfi = 500GeV. 

The SUSY contribution should be compared with the present constraints on the muon 
MDM from experiment and theoretical calculations, which are given in eqs.(|lj) and (0). 
Combining them, we obtain a constraint on the SUSY contribution to the muon MDM, 
Aa^ t USY , which is given by 

- 9.0 x 1(T 9 < AaJ USY < 19.0 x 10" 9 (90% C.L.). (43) 

In Fig. we show the contour of tan/5 which gives the threshold value of the present 
constraint on Aa^ USY given above (i.e. Aa^ USY = —9.0 x 10~ 9 and Aa^ USY = 19.0 x 10 -9 ). 
Here, we choose m^n so that |A^ t USY — 5.0 x 10 9 1 is minimized (where 5.0 x 10~ 9 is the 



center value of the constraint (f43|)). Thus, Fig. |8] should be regarded as a constraint on 
the [iu~ m G2 plane for a fixed values of m^L and tan {3. Notice that if we assume a larger 
value of tan/?, SUSY contribution exceeds the present limit on the muon MDM in wider 
regions. 

Before closing this section, we point out the fact that the contour in Fig. || is not 
symmetric under [in — > —fin- This is because the center value of the constraint given in 
inequality fl43|) is 5.0 x 10 -9 , which deviates from zero. Therefore, constraint fl43|) prefers 
positive value of Aa^ USY , and hence we have severer constraint for fi H < 0. 

6 Discussion 

In this paper, we have investigated the SUSY contribution to the muon MDM by regarding 
all the parameters in MSSM as free parameters. Especially when tan/? is large, the SUSY 
contribution is enhanced, and some parameter region of MSSM is excluded not to conflict 
with the present constraint on the muon MDM. Furthermore, even in the case where tan/? 
is not so large (tan/?^ 10), Aa^ USY may become comparable to the present limit on the 
muon MDM, if the masses of the superparticles are quite light (see Fig. ^a). 
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In MSSM, large tan (3 scenario is an interesting issue that has been received attentions 
in recent years. One of the motivations of large tan (3 is the unification of the masses of 
bottom and tau in SUSY GUT [2{J. That is, in SUSY GUT where the Yukawa coupling 
constants for bottom, y b , and tau, y T , are unified at the GUT scale, the Yukawa coupling 
constant of bottom (or top) is claimed to be significantly large in order to have the 
observed value of the bottom mass. Thus, for the successful unification of y b and y T , 
large value of tan f3 is preferred. (Another solution is to assume tan (3 ~ 1 so that the 
Yukawa coupling constant for top, y t , becomes large.) SUSY GUT based on SO(10) may 
give us another motivation of large tan/3 |21 |. In a simple SO(10) GUT, all the Yukawa 
coupling constants (especially, y b and y t ) are unified at the GUT scale. In this case, tan f3 
as large as m t /mb ~ 50 is required in order to make the hierarchy between the top and 
bottom masses. Furthermore, in some model in which the masses of the light fermions 



are generated radiatively, we need large value of tan/5 ||22|| . The new Brookhaven E821 
experiment will be a powerful test for such types of large tan (3 scenarios. 

Due to the fact that the SUSY contribution to the muon MDM strongly depends on 
tan/?, we may be able to use the muon MDM for the determination of tan/3, especially 
for the large tan /3 case. That is, by future experiments, in particular by NLC, we will be 
able to measure the masses of the superparticles accurately, and it can hopefully fix most 
of the parameters on which the muon MDM depends. Then, precise measurement of the 
muon MDM will give us an useful information about tan (3. 

Comparison of our results with those based on minimal supergravity |Hj may be 
interesting. In both cases, the SUSY contribution to the muon MDM may become 
0(1CT 8 — 1CU 9 ) if tan/3 is large. However, in our result, we can see several interest- 
ing behaviors which hardly occur in the case of minimal supergravity. That is, if we go 
away from the assumption of the universal scalar mass, a cancellation may occur among 
several diagrams when the mass splitting of left- and right-handed smuon is large. Fur- 
thermore, in the case where the SUSY invariant Higgs mass /j,h is quite larger than the 
SUSY breaking parameters, diagram (Nl) in Fig. [| becomes significant, resulting in the 
enhancement of Aa^ USY . 

Finally, we would like to comment on the case with the flavor mixing in the slepton 
mass matrices. In particular, even in the case of minimal supergravity, the sfermion mass 
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matrices receive renormalization effects from the physics much above the electroweak 



scale, such as the right-handed neutrino multiplets Jig, [L3| or GUT |L9]|, resulting in non- 
vanishing off-diagonal elements of the slepton mass matrix. If the off-diagonal elements of 
the slepton mass matrices are substantially large, all the sleptons contribute to the muon 
MDM, as we mentioned before. However, for the case where the flavor mixing exists 
only in left- or right-handed lepton mass matrix, the previous arguments are almost 
unchanged. If both left- and right-handed slepton mass matrices have large off-diagonal 
elements, situation changes. Especially, in this case, Yukawa coupling constant of tau 
can contribute to the muon MDM through the Feynman diagram like (Nl) in Fig. [|, and 
hence the muon MDM may be enhanced. 

Detailed analysis of this case is quite complicated since the muon MDM depends on a 
large number of parameters. Thus, we only discuss the case where the diagonal element 
of the left- and right-handed sleptons, rh 2 L and m 2 R , are proportional to unit matrix; 

a = m L ^rh = 171 r (^ : n °t summed). First, we consider the case where one of rh 2 L or 
rh 2 R has off-diagonal element. In this case, the results of the previous analysis are almost 
unaffected. For example, even if 23/^1, 22 = 0-5 ( or 23/^?? 22 = 0-5), the correction 
to Aa^ USY is less than ~ 10%. If both m 2 L and m 2 R have large off-diagonal elements, 
Aa^ USY may receive a large correction. Numerically, when 23/ r ^i 22 ~ ™r 23/ ^r 22 ~ 
0.2, the correction is 0(10%). The correction gets larger as the off-diagonal elements 
increase. 

The new Brookhaven E821 experiment will give a strong impact on SUSY models. 
By the experiment, the muon MDM is expected to be measured with accuracy about 
0.4 x 10 -9 . Furthermore, the uncertainty in the theoretical prediction, which mainly 
comes from hadronic contributions, is hoped to be reduced by several experiments like 
VEPP-2M, DA$NE and so on. On the contrary, we may have the SUSY contribution to 
the muon MDM to be of order O(10 -9 ) even if all the superparticles are heavier than, 
say, 300GeV (see Fig. |]b) in which case we cannot detect the superparticles even by NLC 
with t/s = 500GeV. Therefore, we may be able to have a signal of the superparticles by 
using the muon MDM even if the superparticles are out of the reach of the forthcoming 
high energy colliders. 
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A Functions In and Jn 

In this appendix, we show some useful formulae for the functions Jjy and Jn, which are 
defined as 

Mml---,m%) = J { i k )H{k 2_ mj) k .. {k 2_ m% y (^) 
The signs of the functions In and .In are given by 

(-l) N I N {ml---,m 2 N ) > 0, (46) 
(-l) N+1 J N (ml,...,m 2 N ) > 0, (47) 

The functions In and Jjv-i are related as 

I N {ml,---,m 2 N ) = — T {I N _i(ml,---,m 2 N _ 1 ) - I N -i{ml, ■ ■ ■ , m 2 N )}, (48) 
m\ — m N 

and the explicit form of I2 is given by 

h(m l m l) = __L i^L-m ( T 4) + In (&)) . (49) 

v 1 21 167T 2 \ml - ml \A 2 J ml -ml V A /J 

Notice that the function I2 is logarithmically divergent, and hence I2 defined in eq. (f49[ ) 
depends on a cut-off parameter A. However, In (N > 3) which is iteratively defined by 



using eq.(f48"D is independent of A, as it should be. In addition, Jn is related to In and 
Ijv-i as 

J N (ml, ■ ■ ■ , m 2 N ) = I N -i{m 2 , • - • , m 2 N _{) + m 2 N I N (m\, ■ ■ ■ , m 2 N ). (50) 
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In the case where all the masses mi - mjv are almost degenerate, it is convenient to 
use the Taylar expansion of 1^. Define 

— 2 2 

m — mf , .„ . 

ei= _ 2 4 * = 1- AO, 51 

with m being an arbitrary mass scale, then ijv is expanded as 

JjvK,-..^^) = 167r2 ^ 2(J v_2) £ (iV + - _ 2 )(iv + p - 1) 

x E 'i ---'x (^>3), (52) 

jiH hjjv=p 

and for N = 2, 

h(<ml) = _^ { ln (g + 1} + ^t^Vnx^ (53) 

Notice that eqs.(^) - (|53|) are useful for numerical calculations. 

Furthermore, the function Jjy has mass dimension (4 — 2N). Therefore, we obtain 

±{\ 2 - N I N (\ml---,\m 2 N )} = 0, (54) 

which reduces to 

N 

(2 - N)I N (ml, • • • , m%) + ^ rn 2 I N+1 (ml, • • • , mf , mf , • • • , m^) = 0. (55) 

i=i 

Similar formula can be obtained for 

N 

(3 - N)J N (ml, ■ ■ -,m 2 N ) + ^m 2 J N+ i(m\, ■ ■ ■ ,m 2 ,m 2 , ■ ■ - ,m 2 N ) = 0. (56) 

i=i 
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Figure caption 



Figure 1: Feynman diagrams which give rise to the muon MDM in the mass eigenstate 
basis. The external lines represent the muon (straight) and the photon (wavy). 

Figure 2: Feynman diagrams which give rise to the muon MDM in the mass insertion 
method. 

Figure 3: The SUSY contribution to the muon MDM, Aa^ USY , in the iiH~ m G2 plane. The 
right-handed smuon mass is taken to be m^R = lOOGeV. We take tan/? = 30, and the 
left-handed smuon mass is (a) lOOGeV, (b) 300GeV and (c) 500GeV. The numbers 
given in the figures represent the value of Aa^ USY in units of 1CT 9 . 

Figure 4: Same as Fig. ^| except for m^ R = 300GeV. 

Figure 5: Same as Fig. |3| except for m^ R = ITeV. 

Figure 6: The SUSY contribution to the muon MDM, Aa^ USY , in the /x#-mG2 plane. The 
right-handed smuon mass is determined so that Aa^ takes its minimal value. We 
take tan/3 = 30, and the left-handed smuon mass rrijXL is taken to be (a) lOOGeV, (b) 
200GeV, (c) 300GeV and (d) 500GeV. The numbers given in the figures represent the 
value of Aa„ Y in units of 10~ 9 . 

Figure 7: Constraints on the /i^-mQ plane for tan (3 = 30 from the negative searches for 
the neutralinos and the charginos. The numbers on the figure represent the lowerbound 
on the chargino mass in units of GeV. The contour with mo2 < 45GeV corresponds to 
the constraint from LEP [16, 17 1. 



Figure 8: Contours which gives the threshold value, i.e. Ac^ USY = -9.0 x 10" 9 (dotted 
line) and Aa^ USY = 19.0 x 10~ 9 (solid line). The right-handed smuon mass m^R is 
determined so that |Aa^ USY — 5.0 x 10 9 1 is minimized. The values shown in the figures 
represent those of tan /3, and we take the left-handed smuon mass rrijxL to be (a) lOOGeV, 
(b) 200GeV and (c) 300GeV. 
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